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The Schwinger model is studied in a finite lattice by means of the P-representation. The vacuum energy, mass 
gap and chiral condensate are evaluated showing good agreement with the expected values in the continuum limit. 



1. Introduction 

A useful formulation of gauge theories, both 
from the conceptual and methodological point of 
view, is the one in terms of gauge invariant ex- 
citations or string-like objects. The so-called P- 
representation 0], consisting of a Hilbert space 
of path labeled states, has been used on the lat- 
tice to perform analytical Hamiltonian calcula- 
tions. A cluster approximation allowed to pro- 
vide qualitatively good results for the (2 + 1) 
QED @ and the (3 + 1) QED § with stag- 
gered fermions. A description in terms of paths or 
strings, besides the general advantage of only in- 
volving gauge invariant excitations, is appealing 
because all the gauge invariant operators have a 
simple geometrical meaning when realized in the 
path space. However, the computational method 
implemented, up to now, on a formally infinite 
lattice, has the serious drawback of the explo- 
sive proliferation of clusters with the order of 
the approximation. In order to tackle this dif- 
ficulty we propose in this paper to explore the 
previous method implemented now on a finite lat- 
tice. As a first test, we choose the simplest lat- 
tice gauge theory with dynamical fermions, the 
Schwinger model or (1+1) QED. This massless 
model can be exactly solved in the continuum and 
it is rich enough to share relevant features with 
4-dimensional QCD as confinement or chiral sym- 
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metry breaking with an axial anomaly (J]. For 
this reason it has been extensively used as a lab- 
oratory to analyze the previous phenomena. The 
lattice Schwinger model also become a popular 
benchmark to test different techniques to handle 
dynamical fermions 

This article is organized into four sections. In 
section 2 we show the formulation of the model in 
the P-representation. The electric and interaction 
components of the Hamiltonian operator are real- 
ized in this basis of "electromeson" states. In sec- 
tion 3, first, we describe the finite lattice Hamilto- 
nian approach. Second, we show the calculation 
of the ground-state energy, the mass gap and chi- 
ral condensate. These results are discussed in the 
concluding section. 



2. Schwinger Model in the Lattice P- 
Representation 

The P-representation offers a gauge invariant 
description of physical states in terms of kets 
| P >, where P labels a set of connected paths 
P$ with ends x and y in a lattice of spacing a. 
In the continuum, the connection between the 
P-representation and the ordinary representation 
("configuration" representation), in terms of the 
fcrmion fields ip and the gauge fields U^(x) = 
exp[iea^4 M (a:)], can be performed considering the 
natural gauge invariant object constructed from 
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them: 



(i) 



where U(P%) = exp[ie Y\ (eP A?] ( I = (x,^) de- 
note the links ). 

The immediate problem we face is that $ is not 
purely an object belonging to the "configuration" 
basis because it includes the canonical conjugate 
momentum of i/j, tp'. The lattice offers a solution 
to this problem consisting in the decomposition of 
the fermionic degrees of freedom. Let us consider 
the Hilbert space of kets | ipt > "^d, A^ >, where u 
corresponds to the up part of the Dirac spinor and 
d to the down part. Those kets are well defined 
in terms of "configuration" variables (the canon- 
ical conjugate momenta of ipd and t/>t are ip d and 
ip u respectively.) Then, the internal product of 
one of such kets with one of the path dependent 
representation (characterized by a lattice path Pf 
with ends x and y) is given by 



$(py)=<P™\ rPl^A, 
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^UxMpy^iy), 



(2) 



where i and j denote a component of the spinor u 
and d respectively. Thus, it seems that the choice 
of staggered fermions is the natural one in order 
to build the lattice P-representation. Therefore, 
the lattice paths -Pjf start in sites x of a given 
parity and end in sites y with opposite parity. 
The one spinor component at each site can be 
described in terms of the Susskind's x( x ) single 
Grassmann fields . The path creation operator 
$q in the space of kets {| P >} of a path with 
ends x and y is defined as 



^Q=xHx)U{Ql)x{y). 



(3) 



Its adjoint operator <tg acts in two possible ways 
Q: annhilating the path Q v x or joining two exist- 
ing paths in | P > one ending at x and the other 
starting at y. 

The Schwinger Hamiltonian is given by 



H 



-W 



W = W E + AW} , A = 



(4) 



W z = + H) (*< = x\x)U n (x)x{x + n)) 

e 

where x labels sites, £ = (x, n) the spatial links 
pointing along the spatial unit vector n, Eg is the 
electric field operator, the kets | P > are eigen- 
vectors of this operator 



E t \P >=N e (P)\P >, 



(5) 



where the eigenvalue Ni(P) is the number of 
times that the link I appears in the set of paths P. 
The are "displacement" operators correspond- 
ing to the quantity defined in (0) for the case of 
a one-link path i.e. P = £. The realization of 
both Hamiltonian terms in this representation is 
as follows jjj: 

By (|^) the action of the electric Hamiltonian is 
given by 



W E \P>=Y,Nf(P)\P> 



(6) 



The interaction term Wi is realized in the loop 
space as 



Wi | P >= y £ d e{P,l) | P-£> 



(7) 



where the factor e(P, t) is or ±1 dictated by the 
algebra of the operators. The different actions of 
operators $f over path-states \P(i) > and their 
corresponding e(P,£) are schematically summa- 
rized in FIG.l. 



3. Finite Lattice Hamiltonian Method and 
Results 

Our method of calculation works assuming a 
lattice of some fixed even number of sites N 
and periodic boundary conditions (PBC). Start- 
ing with the zero-path state | > (infinite cou- 
pling vacuum), then a collection of new states 
| Pi > are generated by applying successively the 
non-diagonal Wi interaction Hamiltonian opera- 
tor - whose action is to add or to eliminate links 
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Figure 1. A summary of the different actions of 
operators i from x to y, applied over path- 
states | P > with their respective e. The original 
path is on the top. The resulting paths \P' > 
are plotted below: (a) represent the addition of a 
disconnected link, (b) the union of 2 disconnected 
pieces, (c) the separation of a connected piece into 
2 pieces and (d) the annihilation of a one-link 
path. 

to to the path P,_i as it was described in the pre- 
vious section - up to order K. The traslational 
symmetry can be exploited in order to reduce the 
dimension of the space tacking only one represen- 
tative P a of each class of translationally equiva- 
lent paths {Pi}. The Hamiltonian matrix, with 
all the transitions between the different states 
\P a >, is then built for the scalar and vector 
sectors and their eigenvalues u>i are numerically 
evaluated. 

In order to perform the generation and recog- 
nition of diagrams (the elementary lattice paths) 
as well as the computation of transitions be- 
tween them, we resorted to the PROLOG lan- 
guage which is very suitable to carry out the sym- 
bolic manipulations. 

The calculations of the ground-state energy, 
mass gap and chiral order parameter were per- 
formed on lattices ranging from size N = 2 to 
N = 16 and at least up to order K = TV in each 
case. Our aim is to extrapolate these results to 
the continuum limit: N — > oo, a — > (A — ► oo.) 

It is clear from the plots (FIGS. 2 to 5) that the 
lattice results show convergence to the expected 
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Figure 2. The ground-state energy density over 
2A 1 / 2 for orders K = 1, 2, . . . , 14 on a lattice of 
size N = 14. 



continuum values. This convergence is, however, 
non-uniform and for A large enough the plots 
show deviation from the continuum values al- 
though the region of assimptotic regime becomes 
larger when the size is increased. It is patent 
that for a fixed lattice size N the best results for 
the vacuum energy and the chiral condensate are 
achieved for order N — 2. This appears to be 
the order at which the finite size effects are min- 
imized. This is not the case with the mass gap 
which always gets closer to the continumm value 
when the order increases. 

GROUND STATE ENERGY 

In the continuum limit the ground-state energy 
density is known exactly [gj : 

lim = — = -0.3183. (8) 

When the order increases loq/(2N\) tends to a 
fixed value. For a fixed size N the closer value to 
(||) is given by order N—2. The value for size N = 
16 and order K = 14 at A = 1000 is lo /(2NX) = 
0.31844, so the discrepancy from the exact value 
is less than 0.05 %. The approximations converge 
with considerable rapidity. FIG. 2 shows for 
orders K — 1, 2, . . . , 14 for A ranging from to 
100 on a lattice of size N = 14. 

In order to obtain a result in a consistent way 
we compute the energy for two large values of 
A, for three correlative large orders and for three 
correlative large sizes. Then, for fixed size and 
order we first extrapolate to A = oo assuming 
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the behaviour a + b/X. Second, for fixed size we 
extrapolate to infinite order assuming exponential 
dependence. Finally we extrapolate to infinite 
size assuming exponential behaviour. The results 
are given in TABLE 1. The error using lattice 
sizes up to N = 16 is 0.17% . 

MASS GAP 

The mass gap for the massless continuum 
Schwinger model can be computed exactly M : 



M c 



1 



= 0.564, 



(9) 



(10) 



e c TrVa 

The lattice mass gap is computed as: 
M u>i — Ll>q 
e ~~ 2VA 

Comparing our results with those of Crewther 
and Hamer jjj obtained by a similar method, 
although they use a different representation 
( Jordan- Wigner transformation), we find com- 
plete agreement for given values of N and K . 
When we reach larger N we observe that the value 
of the mass gap improves substantially. For in- 
stance, in FIG. 3 we show a plot of the mass gap 
for TV = 10 for several orders. As it can be seen 
in the region 10 < A < 30 , the mass gap values 
decrease with the order and the size approaching 
the continuum result. Given the non-uniformity 
of the convergence it is more difficult to extrap- 
olate to the limit although values « 0.579 are 
obtained at the modest size of N = 8. 



CHIRAL ORDER PARAMETER 

An interesting quantity to compute is the vac- 
uum expectation of the chiral condensate per- 
lattice-site < XX >, defined as 



xx = ^ r E(- 1 ) x K t ( x )-^ x )]' 



(11) 



where N s is the number of lattice sites. The 
corresponding operator is realized in the P- 
representation and thus we get for the chiral con- 
densate: 



xx\p >= ( 



1 



A, 



\P>, 



(12) 




Figure 3. The mass gap over 2A 1 / 2 vs. A for 
orders K — 1,2,.., 10 on a lattice of size N = 10 . 



where A/p is the number of connected paths in 
P. To compute < XX > the W Hamiltonian is 
modifyed as 



W' = W+^ £(-lHxt( x ),x(x)], 



(13) 



where a is an arbitrary parameter. Thus, < XX > 
is obtained in the standard way as 



ay. 

< XX >= ~Q^\ a =°- 



(14) 



The massless continuum Schwinger model un- 
dergoes a breaking of chiral symmetry with 



< ifitp > je 



27T 3 / 2 



0.15995, 



(15) 



where 7 is the Euler constant. This non-zero 
value of the chiral condensate is one of the main 
efects of the axial anomaly. 

In FIG. 4 we report the value of the chiral 
condensate per-lattice-site for lattice sizes rang- 
ing from N = 2 to N = 10. FIG. 5 shows this 
chiral order parameter for different lattice sizes 
up to order K = N — 2 for each size. 
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Table 1 



Ground state energy for different lattice sizes and orders at A = oo. 









N = 12 


N = 14 


N = 16 


N = oo 


K = 


N - 


4 


-0.317075 


-0.317723 


-0.318023 




K = 


N- 


3 


-0.318010 


-0.318249 


-0.318348 




K = 


N — 


2 
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Figure 4. The chiral condensate per-lattice-site 
times A 1 / 2 vs. A for orders K = 1, 2, . . . , 10 on a 
lattice of size N = 10 . 



Figure 5. The chiral condensate per-lattice-site 
times A 1 / 2 vs. A on lattices of size N = 4,6, 8, 10 
for order K = N — 2 respectively. 



Notice that the results in the weak coupling 
region converge to the corresponding continuum 
value (|l5]) as long as N increases while for a fixed 
N the value improves with the order K till the 
value K = N — 2 is reached. 

4. Conclusions and Final Remarks 

Our general proposal is to to show that the P- 
representation is a valuable and alternative com- 
putational tool for gauge theories with dynam- 
ical fermions. In particular, in this work, we 
wanted to test the Hamiltonian approach on fi- 
nite lattices. With tis aim, we chose the sim- 



plest model: (1+1) QED. This also enables us to 
compare with the corresponding numerical simu- 
lations jllj using the Lagrangian counterpart of 
the P-representation or the socalled worldsheet 
formulation pd| . This comparison shows that , 
for this case of one spatial dimension, the Hamil- 
tonian method is less time consuming. 

The results are very good and confirm the be- 
lief of Hamer et al || in obtaining with consider- 
able accuracy the observables working on lattices 
of moderate size. Consequently, this procedure 
is appealing because one can run all the needed 
computations in small machines obtaining quite 
fair results. 
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Finally, we would like to stress (once more) that 
our aim was not to present another solution to 
the Schwinger model, but, to test an alternative 
general approach to tackle dynamical fermions. 
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